Quantum searching with continuous variables 
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A fast quantum search algorithm for continuous variables is presented. The result is the quan- 
tum continuous variable analog of Grover's algorithm originally proposed for qubits. A continuous 
variable analog of the Hadamard (i.e., Fourier transform) operation is used in conjunction with in- 
version about the average of quantum states to allow the approximate identification of an unknown 
quantum state in a way that gives a square-root speed-up over search algorithms using classical 
continuous variables. Also, we show that this quantum search algorithm is robust for a generalised 
Fourier transformation on continuous variables. 
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Quantum systems can register and process informa- 
tion in ways that classical systems cannot. As a result, 
it is possible for quantum computers to perform certain 
computational tasks faster than any classical computer 
It is becoming increasingly clear that at the heart 
of quantum computation lies two basic quantum phenom- 
ena, one is quantum interference and the other quantum 
entanglement. The real upsurge of interest in this field 
came after Shor's flpfl remarkable discovery of an algo- 
rithm for factoring large numbers |pT| . Subsequently, a 
fast quantum search algorithm has been discovered by 
Grover [Q, which takes 0(\/~N) steps instead of N steps 
to search an unmarked item in a unsorted list of N en- 
tries. Later the optimality of this search algorithm was 
proved |ll|[l4]]. In particular, it was shown that this 
search algorithm can use almost any unitary transfor- 
mation on qubit states Jl5| . Further, it has even been 
argued that this algorithm would work even without en- 
tanglement (though at a cost in resources) |1(J . 

These algorithms are usually implemented on quantum 
systems whose observables have discrete spectra, such as 
a collection of two-level atoms, ions, or spin-i particles 
(called qubits). However, there are other classes of quan- 
tum systems whose observables form continuous spectra. 
So, it is important to know how these algorithms can 
be generalised for continuous quantum variables? With 
the recent advances in our ability to manipulate contin- 
uous quantum information in teleportation fl7j| , error- 
correction codes Jl|,[n| and its feasibility of implementa- 
tion using linear devices |2(i(| , it is natural to ask whether 
one can provide some quantum algorithms that might 
be implemented on a continuous variable quantum com- 
puter. In fact, the usefulness of quantum computation 
over continuous variables has been recently emphasised 
plf . It has been shown that universal quantum compu- 
tation over continuous variable is not only possible, but 
could be effected using simple non-linear operations with 
coupling provided solely by linear operations . These 
operations form a universal set of quantum gates for 



continuous variables allowing 'quantum floating point' 
arithmetic. While discrete quantum computation can be 
thought of as the coherent manipulation of qubits, con- 
tinuous quantum computation can be thought of as the 
manipulation of 'qunats', where the qunat (pronounced 
as 'Q nat') is the unit of continuous quantum informa- 
tion. 

In this letter, we propose a fast quantum search al- 
gorithm with continuous variables. Here a continuous 
variable can be anything, e.g., position, momentum, en- 
ergy (unbounded) or amplitudes of the electromagnetic 
field. With the help of the Fourier transform (viewed as 
an active operation) on a continuous basis state (analo- 
gous to the Hadamard transform in the case of qubits) 
and a suitably generalised inversion operator, we con- 
struct a search operator which can be implemented on 
a continuous variable quantum computer. The inversion 
operator requires the projection operator for continuous 
basis states which we discuss. We show that the ap- 
plication of the compound operator takes 0(^f~N) itera- 
tions to search an unmarked item in a list of N entries. 
Further, we generalise our quantum searching with con- 
tinuous variables to one based on a generalised Fourier 
transformation and which still gives a square root reduc- 
tion in the number of steps. This shows that the quan- 
tum search algorithm with continuous variables is robust 
to the choice of arbitrary Fourier transformations. We 
also discuss the robustness of the search algorithm if one 
uses strongly peaked normalisable states instead of ideal 
infinite-energy position eigenstates. 

Here, we discuss how to perform a quantum search 
algorithm using continuous variables. First we need 
to map a conventional discrete search problem into a 
continuous variable context. Suppose we have a func- 
tion f(k) : K — > {0,1} defined on a domain K with 
k 6 K = {1, 2, . . . , N}. This function has a non-zero 
value equal to 1 for some element k = kf and is for 
all other elements in the set K. Our task is to discover 
the value of kf given the ability to apply the function / 
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to inputs or superpositions of inputs, and given no fur- 
ther information about the function f(k). In order to 
implement this in a quantum computer with continuous 
variables we require a collection of n qunats. The state 
vector of each qunat belongs to a Hilbert space of in- 
finite dimension. Since we have an infinite number of 
basis states, we cannot map each basis state within the 
Hilbert space onto each entry in the set K . (This would 
be a many-to-one mapping.) One could avoid this prob- 
lem by choosing a subspacc of the full Hilbert space with 
N disjoint regions of the spectrum. However, we do not 
discuss this approach in detail here. 

Let us consider a collection of n continuous variables 
whose Hilbert space is spanned by a basis of states 
| a;) = X2, ■ ■ ■ , x n ), satisfying the orthogonality con- 
dition (x\x'} = S(xi — x[) ■ ■ ■ S(x n — x' n ) — 8{x — x'). For 
example, one can consider a compact region of the state 
space divided into TV equal subvolumes, each with mea- 
sure Ax n , one for each member of the set K. Let Xf be 
the centre of the subvolume corresponding to kf . In the 
context of this continuous variable embedding, executing 
the function / corresponds to adjoining an extra state 
to the system, originally in the state |0), and applying 
an operator Uf : |x)|0) — > |x)|l) if x belongs to the re- 
gion corresponding to kf and |x)|0) — > |x)|0) otherwise. 
Clearly, if one samples the region at random by applying 
the operator to a series of random points, it will take 
O(N) calls of the operator to find kf. 

If one exploits the power of quantum superposition and 
entanglement, however, fewer function calls are required. 
Our approach is discussed below. Let us pick an initial 
state in the position basis such as \xi) = \x\, X2, ■ • ■ , x n )i 
for a quantum computer with continuous spectrum at 
random. The final (target) state is given by \xf) = 
\x±, x%, . . . , x n ) f . We need a suitable unitary operator, 
which can take the initial state to the final state. Just as 
we have the Hadamard transformation in discrete com- 
putation, one of the basic operations with continuous 
variables is the Fourier transformation between position 
and momentum variables in phase space. By defining the 
Fourier transformation as an active operation on n qunat 
states \x) we can write it as 

f\x) = -±=Jdye 2i *y\y), (1) 

where xy = Xiyx H h x n y n , \y) = \yi, y 2 , . . . , y n ) and 

both x and y are in the position basis. This has been 
used by one of the present authors (l^j2(| in developing 
an error correction code for continuous variables. This 
Fourier transformation can be straightforwardly applied 
in physical situations. For example, when \x) represents 
quadrature eigenstate of a set of modes of the electromag- 
netic field, J-\x) is simply an eigenstate of the conjugate 
quadrature. 

Suppose, we apply the unitary operator J 7 to a ba- 
sis state \xi), then the relative amplitude of finding the 



system in the target state \xf) is {xf\T\xi) = Tfi = 
e 2txix f 1 _ xherefore, the relative probability of find- 
ing the system in the final qunat states will be given by 
|.F/i| 2 = I/71-". Hence, we have to repeat the experiment 
at least \/\Tfi\ 2 — n n times to successfully obtain the 
state \xf). Here, we prove that search algorithm based 
on continuous variable can take y/n™ steps to reach the 
final state starting from an initial state. (Here, we may 
identify the number of entries N with 7r n .) 

The next operator we need is the unitary operator, 
which can invert the sign of a basis state \x). We can 
define the selective inversion operator for a continuous 
basis \x) as 

I x = l- 2P Ax , (2) 

where P Ax is the projection operator for continuous vari- 
ables. Unlike the discrete case we cannot define the pro- 
jection operator for the basis |x) as P x = |x)(x|, because 
the operator P x is an ill defined and it will not satisfy 
P 2 = P x - The correct projection operator for continuous 
variables is defined |22) as 

i.xq+Ax/2 

Pax = / dx'\x')(x'\ . (3) 

Jx -Ax/2 

The reason for this definition is that we cannot project 
an arbitrary state which is represented in terms of contin- 
uous basis state onto a point to get the exact eigenvalue. 
There will be always a spread within an interval. We can 
only project a state around xq to a selectivity Ax of the 
measuring apparatus. It is not possible to design a device 
to make a perfectly selective measurement of a continu- 
ous variable. The interval [xi,X2] cannot be narrowed 
down, because it will always contains an infinite num- 
ber of eigenvalues p2| . Thus, if we have a wave packet 
the effect of projection is to truncate it around xq within 
an interval Ax. This operator satisfies P\ x = Pax and 
Pax\x) = |x) as expected. With the help of the above 
inversion operator we can construct a compound search 
operator C defined as 

c 1, ./• /, :/•. (4) 

It may be remarked that the selective inversion of the 
target state |x/) can be achieved by attaching an ancilla 
qunat and considering the quantum XOR circuit for con- 
tinuous variables fisj| . If a quantum circuit exists that 
transforms |x)|a) — ► |x)|/(x) + a), then by choosing the 
ancilla state \a) = T\nj r 2) = J dy e tny \y) / 'a/tt™ we can 
selectively invert the state |x) for which f(x) — 1, i.e., 
\x)T\7r/2)^-\x)T\7r/2). 

Let us define a state |x/) = T^\xf). We can show that 
the operator C can preserve the two-dimensional subspace 
spanned by the states \xi) and \xf). First, we show the 
action of C on \xi). This can be expressed as 

C\xi) = \x t ) - AP Axt T ] P Axs T\ Xl ) + 2pP Axf F\ Xl ) , (5) 
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where Pax, = dx\ \x')u(x'\, (x n = x -Axi/2, x l2 = 
xq + Axi/2 ) and likewise for Paxj ■ Using these facts we 
can simplify the above equation to 



C\xi) = (1 -)\xi 

7T" 



Xf2 



dx'f e 



(6) 



Similarly, we can evaluate the action of C on \xf). It is 
given by 



C\ 



(7) 



Thus, the operator C creates superpositions of two qunat 
states just as Grover's operator creates superpositions 
of two qubit states. Once we understand the action of 
C on qunats we can obtain the total number of steps 
required in reaching the target state. Here, we use ge- 
ometric structures from the projective Hilbert space of 
a quantum system to obtain the number of steps in the 
quantum searching. The projective Hilbert space admits 
a natural measure of distance called Fubini-Study dis- 
tance [|3| . This measures the shortest distance between 
any two (not necessarily normalized) states and l^) 
whose projections on V are n( , 1 ) and n(i/>2)) respec- 
tively. This can be defined as 



d 2 (|Vi),IV>2» = 4 1 
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Here, the vectors l^i) and \ip2) can be quantum states 
over continuous variables or over discrete variables. For 
unnormalisable states the above definition still works pro- 
vided it is understood that the norm of the states can be 
made finite. In dealing with position eigenstates we can 
imagine that either the particle is moving in a finite space 
PU so that position eigenstates do not diverge or one can 
use normalisable states having strong peaks around some 
value of the position axis. 

During the quantum searching with continuous vari- 
ables we want to reach a state \xt) from an initial state 
\xj). This means we have to travel a shortest distance 
between these states which is given by d 2 (\xi), \xf)) = 
4(1 — l/w n ). One application of the operator C creates a 
state \xi)^' = C\xi). We calculate the shortest distance 
between the resulting state \xi)^ and the initial state 
\xi). We note that the overlap of these states is given 
by (xi\C\xi) = (1 - A/TT n ){x l \x l ) + 2Ax f /ir n . For large 
database search N — 7r™ is very large and if we assume 
that the measuring device has a narrow selectivity, then 
Axf is also small. Hence, we can neglect the second term 
in the overlap (as it is a product of two small terms). Also 
note that term (xi\xi) is not normalised but nevertheless 
it cancels out in during calculation. With this idea in 
mind we can evaluate the shortest distance between these 



states which is given by d 2 (\xi 



)) « 32/71-™. Thus in 



one application of the search operator C we can move the 



initial basis a shortest distance 0{l/\fK™). Therefore, to 
travel the full distance on the quantum state space we 
need N s number of steps given by 
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This shows that a quantum computer based on qunats 
can take 0(\/n™) applications of C to reach the target 
state which otherwise would have taken 0(jr n ) number 
of steps by the application of T on \x%). Because the 
state is moving along a geodesic each application of C 
rotates the initial state in the right direction. This is our 
quantum search algorithm with continuous variables. 

Instead of position eigenkets one can use strongly 
peaked normalisable state such as 

dx exp[ Jj^fl ]]x) , (io) 
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When e — > 0, the state \rt) becomes a position eigenstate 
\xi). Our algorithm can be practically implemented with 
such states. One can see that the action of the search 
operator C on |r») gives 

4„ , 2 1 
n' 



l 



dx 



)\n) + 

Xf2 



dx'je 



7T™ (27re) n / 4 

(x-x,) 2 /4e 2 



- 2ixx 'fJ*\x' f ) 



(11) 



From the above formula one can see that a single appli- 
cation of search operator moves the initial state |r») a 
distance given by d 2 (\ ri ), |rf } )) w 48/tt". Now, if one 
defines the target state as 



\r f ) 
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dxeM- (X J /)2 ]N) 



(12) 



then one can check that the total (shortest) distance be- 
tween the initial state |rj) and the desired state = 
J^\r f ) is 4[1 - 40(e 2 )/7r™]. Hence, by using (9) the to- 
tal number of steps required to reach the target state is 
N s = O(V^). 

Now we show that the quantum search algorithm over 
continuous variables is robust to some extent. Instead of 
the Fourier transform T if we replace it by a generalised 
Fourier transform (GFT) in the search operator C, still 
the algorithm works, i.e., we do get a square root re- 
duction in the number of steps. We define a generalised 
Fourier transform as an active operation in the position 
basis \x) as 



n/2 



7rsm( 



x dy exp 



sin 9 



[{x 2 +y 2 )cos0- 2xy] 



\V) ■ (13) 



The GFT with a flexible angle 9 gives a physical change of 
the basis \x) by any desired amount p3j. Here, it should 
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be mentioned that 9 > arcsin(l/7r), since for smaller val- 
ues of 9 the assumption sin™ 9 ^> l/w n does not hold. 
The GFT for 9 = 27rm, m being an integer, corresponds 
to no change of basis. The GFT for 9 = n/2 corresponds 
to the Fourier transform defined in (1) (up to a constant 
phase shift equal to nn/A, n being the number of qunats). 
If we apply GFT to an initial basis \xi) then by probabil- 
ity rules of quantum theory we have to perform at least 
0{(n sin9) n ] number of trials to reach a target state \xf). 
We will prove that the generalised search operator acting 
on continuous variables will take Ok/ (tv s'm9) n ) steps to 
reach the final state. 

The search operator with this GFT takes the form 

C {6) =-I Xi F {e)] I Xf F {9) . (14) 

We can see that the action of the generalised search op- 
erator on the initial state \xi) is given by 

4 / i n 

C W \*<) = I 1 - f g1 - avJ to) + 2 \ 1 • ( 15 ) 
(7rsm6')™ y yTT sm 9j 

f Xf2 i 
x / dx' f exp{-^— [(x 2 + x'f) cos6> - 2x i Xf})j : ^\x'f) . 

Jxf! SVO.U 

Similarly, the action of the generalised search operator 
on |x)/ can be calculated. It is given by 

cm ^' = ^'- 2 {Mw (16) 

x J dx' l exp{ — [(x ,2 +x 2 f )cos9~2x' t .x f }}\x^. 

It can be seen that the generalised search operator cre- 
ates linear superposition of qunat states in the search 
process. Now, we can calculate the Fubini-Study dis- 
tances to know how many steps are needed to reach the 
target state. The shortest distance between the states 
\xi) and \x f ) is d 2 (\x l ), \x f )) = 4[l-l/(7rsin0)"]. Notice 
that single application of the search operator C^ s ' moves 
the initial state by a distance given by d 2 (\xi), = 
(32/7rsin#)™. Therefore, to travel a shortest distance 
d(\xi), \xf)) we need N s w Ok/ (7rsin#)™] number of 
steps. 

Thus, using a generalised Fourier transform we have 
proved that there is a square root reduction in the num- 
ber of steps working with continuous variables. As ex- 
pected for an angle 9 = ir/2 we get back the original 
result with the search operator C. This result is simi- 
lar to the recent result of Grover Jig] , where the search 
algorithm for qubits has been generalised for arbitrary 
unitary transformations. 

In conclusion, we have for the first time provided an 
efficient algorithm such as quantum searching to be im- 
plemented on a quantum computer with continuous vari- 
ables. The key elements in this generalisation are the 



Fourier transformation and inversion operators which 
constitute the search operator for qunats in an infinite 
dimensional Hilbert space. We find that a square root 
speed up is possible with quantum computers based on 
qunats. Also, the continuous variable search is possible 
with almost any Fourier transformation. This may be 
practically implemented for any large data base search 
using linear and non-linear optical devices with the role 
for qunats being played by electromagnetic fields. It may 
well be that for large data base searches it is beneficial 
to use continuous quantum variables. 

AKP and SLB acknowledge financial support from EP- 
SRC. 
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